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We construct a new class of plane-symmetric solutions possessing a cur- 
vature singularity which is null and weak, like the spacetime singularity at 
the Cauchy horizon of spinning (or charged) black holes. We then analyse the 
stability of this singularity using a rigorous non-perturbative method. We find 
that within the framework of (linearly-polarized) plane-symmetric spacetimes 
this type of null weak singularity is locally stable. Generically, the singularity 
is also scalar-curvature. These observations support the new picture of the 
null weak singularity inside spinning (or charged) black holes, which is so far 
established primarily on the perturbative approach. 



I. INTRODUCTION 

The Kerr solution represents the geometry of the general stationary, vacuum, spin- 
ning black hole (BH). This solution has an obvious relevance to reality, because realistic 
astrophysical BHs are generally spinning . It is well known that a Kerr BH has a Cauchy 
horizon (CH) - a null hypersurface which marks the boundary of the domain of depen- 
dence for any initial hypersurface in the external world. A similar CH also exists in the 
Reissner-Nordstrom (RN) geometry, which describes spherically-symmetric charged BHs. 
The presence of a CH inside a spinning or charged BH is disturbing, because the Einstein 
equations do not provide a unique prediction for the extension of the geometry beyond the 
CH. However, it is also well known that the CH of RN and Kerr is unstable; Namely, linear 
perturbations of various types develop singularities at the CH, [R|- || suggesting that in a 
generic situation the smooth CH of Kerr or RN will be replaced by a curvature singularity. 

For many years, the nature and exact location of this singularity were not completely 
clear. The prevailing point of view in the last few decades was that the spacetime singularity 
inside BHs is of the BKL [[J type, i.e. spacelike, oscillatory, and tidally-destructive. In the 
last few years, however, there is steadily growing evidence that a singularity of a completely 
different character forms at the CH of spinning [|T~0|JTTf| and charged |12|]- [15|] BHs. This sin- 
gularity is null [}TD| , P~2| , [1~5| , P~5| ] rather than spacelike, and weak |TU| , |T3[ (in Tepler's terminology 
|17|1), rather than tidally-destructive. 



So far, the most direct evidence for the formation of a null weak singularity inside generic 
spinning BHs stems from the nonlinear perturbation analysis of the interior Kerr geometry 
101)111 . It would obviously be important to confirm the perturbative results from alterna- 



tive, non-perturbative, directions of research. Motivated by this, Yurtsever |L8| suggested 
that plane-symmetric spacetimes could serve as an excellent test-bed for further exploring 
and testing the new picture of the BHs' null weak singularity, emerging from the perturba- 
tive approach. Yurtsever based his suggestion on the following argument: If indeed a null 
curvature singularity exists at the CH of generic spinning (or charged) black holes, there 



should exist corresponding plane-wave solutions fT9||20|j which admit a (locally) similar type 



of null singularity. Therefore, understanding the rule of such null singularities in plane-wave 
and colliding plane-waves (CPW) solutions [0,0] may provide important insight into the 
issue of stability of the null CH singularity inside spinning (or charged) BHs. 

In this paper we shall attempt to undertake this goal. We shall first construct an exact 
linearly-polarized (LP) ingoing plane-wave solution which admits a weak curvature singu- 
larity on a null hypersurface. Then we shall analyse the stability of this type of singularity, 
within the framework of (LP) plane-symmetric spacetimes, by a rigorous non-perturbative 
method: First we shall demonstrate the stability with respect to generic ingoing plane-wave 
perturbations (which preserve the plane- wave character of the solution). Then we shall 
introduce outgoing plane-symmetric perturbations as well, and analyse their effect on the 
structure of the singularity. The outgoing perturbations convert the geometry into a (LP) 
CPW spacetime. We shall show that the singularity remains null and weak, though curva- 
ture scalars, which were strictly zero in the original plane-wave solution, generically blow 
up when the outgoing perturbations turn on. (This situation is fully analogous to what is 
known about the black holes' CH singularity: In a spherical charged BH, if the radiation is 
purely ingoing the CH singularity is non-scalar ||16|| , but when one adds outgoing radiation, 
the singularity becomes scalar-curvature - the mass- inflation singularity |12" |. In the case of 
a spinning vacuum BH, the CH singularity is scalar-curvature. [ ]TU| , [2"5] ) We conclude that the 
null weak singularity is locally stable within both frameworks of plane-wave solutions and 
CPW solutions (though the non-scalar character of the null singularity in plane-wave solu- 
tions is unstable within the framework of CPW solutions). This provides a strong support 
to the above-mentioned new picture of the BHs' singularity. 

In Ref. |]18[ Yurtsever presented a certain limiting process which maps a (local neighbor- 
hood of a) generic null singularity into a plane-wave null singularity. Based on this limiting 
process, Yurtsever correctly pointed out that if indeed a generic null singularity exists inside 
black holes, it should be asymptotically similar to a generic null singularity of a plane-wave 
solution. Yurtsever further argued that this plane-wave null singularity should coincide with 
the "singular Cauchy horizon" of the plane-wave solution. But the Killing Cauchy horizons 
of plane-wave solutions are known to be non-generic and unstable within the context of 
global initial-value problem for CPW spacetimes. JEJ This led Yurtsever to the conclusion 
that a null singularity (e.g. inside black holes) must be locally unstable and hence unrealis- 
tic. This conclusion, we argue, is incorrect; Its derivation requires one to make two closely 
related assumptions: 

Assumption 1: The scenario of a formation of a black hole, with a null curvature singularity 
inside it, from regular initial data, can be approximated in a global sense by some plane- wave 
(or CPW) solution with a null singularity, which evolves from regular initial data. 
Assumption 2: Correspondingly, the latter null singularity should be located at the Killing 
Cauchy horizon of the plane-wave (or CPW) solution. 

Both assumptions result from the confusion of local and global aspects of the problem. The 
statement that the black-hole's null singularity is well approximated by a plane-wave null 
singularity holds only on a local basis. This is obvious from the nature of the limiting process 
described in Ref. fllSfl , which assumes a null geodesic 7 that intersects the null singularity: 



It is only the "tubular" immediate neighborhood of 7 (near its intersection point with the 
singularity) which admits the approximate similarity to a plane-wave solutions. (See Fig. 1, 
in which this neighborhood is denoted by N.) Since the initial hypersurface of the black-hole 
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spacetime is remote from the null singularity inside the BH, there is no way to extend TV 
up to this initial hypersurface. Therefore, despite the local similarity of the neighborhoods 
of the null singularities in the two spacetimes, their global features are totally different. In 
particular, there is no reason to relate the BH's null singularity to the plane-wave Killing 
Cauchy horizon. Moreover, when considered as a simplified approximate model for the BH's 
null singularity, the issue of stability of the plane-wave null singularity must be considered 
from the local point of view (see below), and not from the global one based on Cauchy 
evolution from regular initial data. 




FIG. 1. The local initial-value setup for the plane-symmetric spacetime, and its relation to 
the inner structure of the black hole. EH and CH are the event horizon and the singular Cauchy 
horizon of the black hole, respectively. The two characteristic initial null hyper surf aces, u = uq 
and v = vq, are located in the neighborhood N of the intersection point of the CH singularity and 
the outgoing null geodesic 7. 

Our approach to the issue of local stability is based on a local (characteristic) initial- 
value setup. Namely, we introduce a pair of characteristic initial null hypersurfaces (denoted 
Mo and Vq in Fig. 1), both located inside the small local neighborhood iV of a point on the 
null singularity. In such an initial-value setup, the generic CPW solution depends on two 
freely-specifiable functions, one of the ingoing null coordinate v and one of the outgoing 
null coordinate u. (In the special case of a plane-wave solution, one of these functions 
degenerates to a constant.) We shall first construct a specific ingoing plane- wave solution 
which admits a null curvature singularity at v = v s . Then we shall add perturbations to both 
the above mentioned initial functions. These perturbations are generic, though bounded 
(and restricted in some weak sense). We shall show, by means of a non- perturbative 
analysis, that the null singularity at v — v s survives these perturbations. That is, within 
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the framework of CPW spacetimes, the null singularity is locally stable. 

In this local setup, the outgoing initial hypersurface w necessarily intersects the null 
singularity (obviously, if we chop this hypersurface before v = v s , the null singularity will 
not be included in the domain of dependence). In other words, in order to recover the 
null singularity, we must start with initial data at u = u$ which are themselves singular at 
a certain point (u = uq,v = v s ). |24| (This singular point will then be the "seed" of the 
null singularity, which will extend along the hypersurface v = v s .) This disadvantage is 
unavoidable, because, as was explained above, the global plane-wave initial-value problem 
(starting from regular initial data) has little to do with the issue of null singularities inside 
black holes. In physical terms, what we show here is that a singularity (of a certain type) at 
the initial data generically propagates along the characteristic line and thereby produces a 
null singularity. This observation is not trivial: One might conceive that the nonlinearity of 
the field equations will generically prevent the formation of a null singularity, and a spacelike 
singularity will form instead (this is essentially the possibility suggested by Yurtsever [|T%|j). 
Our local stability analysis shows that this is not the case: The null singularity of plane 
waves is locally stable within the frameworks of plane-wave and CPW spacetimes. 

Recently, Ori and Flanagan (OF) p5| used another construction to demonstrate the lo- 
cal genericity of null weak singularities. The construction by OF is more powerful than the 
one presented here, as it demonstrates local (functional) genericity within the framework 
of the fully-generic class of vacuum solutions, without any symmetry. The present work 
is advantageous in certain respects, however: First, the mathematical construction by OF 
heavily depends on analyticity: It is restricted to analytic initial-value functions, and is pri- 
marily based on the Cauchy-Kowalewski theorem. The present construction does not make 
any assumptions about the analyticity of the initial functions (we only demand smooth- 
ness), and is based on the hyperbolic nature of the field equations. Second, the type of 
local asymptotic behavior covered by the present analysis is slightly more general than that 
of OF, as explained in section M. In addition, the generalization to null singularities with 
more realistic types of local asymptotic behavior will be much easier to implement within 
the present framework of CPW spacetimes (this is again a consequence of the lack of any 
reference to analyticity in the present method; see section |V[) . Third, the simplicity of the 
CPW solutions makes the present construction more useful for investigating various features 
of the null singularity. We should also mention another attempt to use a local construction 
for investigating the BH's null CH singularity, made earlier by Brady and Chambers |26| . 
Their construction, however, only addressed the constraint equations on two null hypersur- 
faces, and the compatibility with the six evolution equations was not considered there. In 
particular, the analysis in Ref. [ZB| does not rule out the possibility that a spacelike singu- 
larity will form immediately at the singular "point" [analogous to (uq,v s )] at the future edge 
of the outgoing initial null hypersurface. The present construction (like that of OF) takes 
care of the full set of vacuum Einstein equations. 

Throughout this paper we restrict attention to linearly-polarized (LP) plane-wave and 
CPW solutions, for simplicity. We do not expect the qualitative results to be different in the 
more general case of arbitrarily polarized solutions. The analysis in the arbitrarily polarized 
case is more complicated, but still appears to be feasible. 

This paper is organized as follows. In section [D] we shall construct an explicit ingoing 
LP plane-wave solution with a non-scalar parallelly-propagated (PP) null weak singularity, 
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to which we shall refer as the basic solution. Then, in section [1(1] we shall perturb this 
basic solution by generic (though bounded) LP ingoing perturbations, and show that the 
singularity remains non-scalar, PP, null, and weak. This demonstrates that the non-scalar 
null weak singularity is a generic feature of (LP) plane-wave spacetimes. In section [V] we 
shall add generic (LP) perturbations in the outgoing direction. The geometry is now de- 
scribed by the CPW solution. We shall show that the singularity remains null and weak, 
though generically the outgoing perturbations convert the non-scalar PP curvature singu- 
larity into a scalar-curvature singularity. Finally, in section [V| we shall discuss the extent 
and implications of our results. 

II. BASIC SINGULAR PLANE- WAVE SOLUTION 

The LP plane-wave spacetimes can be described by the line element 

ds 2 = -du dv + e~ u ^ [e v ^dx 2 + e~ v ^dy 2 ] . (1) 

The only non-trivial vacuum Einstein equation is 

2U, W - U v 2 -V v 2 = 0. (2) 

Since there are two non-trivial functions [U(v) and V(v)] and one constraint, Eq. (2), plane- 
wave solutions are described by one arbitrary function of v. To simplify the calculations, 
we take the freely-specified function to be U(v). Equation (2) then becomes an integral for 
V(v). [When the freely-specified function is taken to be V(v), Eq. (2) becomes a non-linear 
differential equation for U(v).] 

We wish to construct an explicit solution in the range vq < v < v s which develops a 
null weak curvature singularity at v = v s . In general, such a solution will be obtained from 
any function V(v) which is smooth at v < v < v s and continuous at v = v s , but with 
V v diverging at v — v s . We shall now construct a simple explicit solution of this type. To 
shorten the notation, we take v s = and define 



w = \ v \V n = (- v yi n , (3) 

where n is any constant > 2 . (Note that with the choice v s = 0, both v and Vq are negative.) 
We now take 

U(v) = - (n - 2) a 2 w 2 = U (v) , (4) 
where a is a positive dimensional constant. Solving Eq. (2), we obtain 



V(v) = (n - 2) 



1/2 



uiu (l — a 2 w 2 ^j + arcsin (aw) =Vq(v). (5) 



(For concreteness we took here the positive root for V v , and set the integration constant 
such that V vanishes at v = 0. This causes no loss of generality: Adding a constant to V 
or changing its sign amounts to rescaling the coordinates x and y or interchanging them, 
respectively.) In order to ensure the validity and regularity of V(v) in the entire range 
^ v < 0, we take < a < |fo| _1 ^ n , so that < aw < 1. 
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We shall refer to the explicit plane-wave solution (4,5) as the basic solution. In terms 
of our characteristic initial-value setup, this plane wave solution may be viewed as evolving 
from the initial data (4,5) at the characteristic hypersurface u = Uo, with trivial (i.e. u- 
independent) initial data along the other hypersurface, v = v . 

An expansion of Eq. (5) near v = yields 

V = 2(n - 2)aw + O (w 3 ) . (6) 

For any function F(w), we find 

F v = _ n -i\ v \i/n-i lw } F vv = n- 2 [\v\ 2 / n - 2 l ww -(n-l)\v\^ n ~ 2 l w ] . (7) 

Therefore, at v — the derivatives of V and U with respect to v diverge: 

V v ^- [2n- 1 (n-2)a] \v\ l/n ~ l , U v oc v 2 ^ 1 . (8) 

(This divergence has an invariant meaning, because v is the affine parameter for the null 
geodesies of constant x, y, u .) The second-order derivatives diverge even faster: 



V. 



vv — 



2n- 2 {n- l)(n-2)ol \v\ 1/n ~ 2 , U vv oc v 2/n ~ 2 . (9) 



As a consequence, various components of the Riemann tensor diverge at v — 0. For example, 
one finds that [] | 

Rvxvx = \e V ~ U {-V vv + U V V V ) , (10) 

and thus near v — 0, 

R vxvx = ~\v vv = n- 2 (n - l)(n - 2)a\v\ 1 / n - 2 . 

Since all curvature scalars vanish in plane-wave spacetimes, there is no scalar-curvature 
singularity at v = 0. One can easily verify, however, that parallelly-propagated (PP) com- 
ponents of the Riemann tensor diverge along null and timelike geodesies intersecting v = 0. 
For example, along a null geodesic of constant convenient PP tetrad is 

e a = e (U-V)/2 5 a > g a = e (U + V)/2 g« > £ = % , e « = ^2 ^ . (11) 

The PP tetrad component 

e'VeVi? a a = 2e {u - y) R 

°D c i c i - 1 L ap"/o - 1 l ra»a; 



This expression is obtained from a direct calculation of the Weyl tensor, which is equal to the 
Riemann tensor in the vacuum spacetimes considered here. 

2 The most divergent components of Riemann are R VX vx and R vyvy (R V xvy vanishes identically). 
The expression for R vyvy is the same as that of R VX vx, except that V is replaced by —V. (The 
same relation holds for the corresponding PP components discussed below.) 
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diverges like \v \ 1 / n ~ 2 . 

Thus, there is a PP curvature singularity at v — 0. However, because U — V — there, 
the line element (1) is well-define (i.e. finite and non-degenerate) even at v = 0. It then 
follows that the curvature singularity at v = is weak (in Tiplers terminology [jT7||). We 
conclude that the hypersurface v = is a null, weak, PP curvature singularity. Note also 
that no other singularity occurs in the range v o < v < , because U and V are smooth 
(C°°) functions of w. 



III. STABILITY TO INGOING PLANE- WAVE PERTURBATIONS 

In this section we shall analyse the stability of the basic plane-wave solution (4,5) within 
the framework of (LP) plane waves. The basic solution is characterized by the initial func- 
tions Uq{v) and V (v), Eqs. (4,5). We shall now add a small (though finite) perturbation 
AU(v) to the initial function Uq(v): 

U(v) = U (v)+AU(v) . (12) 

In terms of the characteristic initial- value setup, this amounts to perturbing the initial data 
at the characteristic hypersurface u = uo, while leaving the trivial (u- independent) initial 
data at v = v o unchanged (apart from a trivial constant shift, to allow for continuity at the 
intersection point of the two initial hypersurfaces) . We assume that the perturbation AU(v) 
is a smooth function of v in the range v$ < v < 0. By virtue of the constraint equation (2), 
this perturbation of U(v) will lead to a corresponding perturbation in V(v ), which we now 
analyse. We first write Eq. (2) as 

V 2 = 2U m - U/ = V / + [2AU m - AU/ - 2U , V AU V 

Converting the derivatives of V, Vo, and Uo (but not of AU) from v to w, recalling v iW = 
—nw 11 ' 1 , and selecting the positive root (as before), we obtain an expression of the form 

V w = 2a(n - 2) [l - a 2 w 2 + hAU v w n + (b 2 AU/ + b 3 AU tVV )w 2n - 2 ] ^ , (13) 

where bi, bi and 63 are constants that depend on a and n. In the basic solution (AU = 0), 
the term in brackets is strictly positive in v < v < (recall that \aw\ < 1). Therefore, 
bounds Bi,B 2 (which may depend on a, n, and Vq) exist such that the term in brackets will 
be strictly positive in v < v < for any AU (v) satisfying 

|AEg<-B! , \AU m \<B 2 (14) 

throughout this range. For all perturbations satisfying the inequalities (14), V w is well- 
defined and smooth in Vq < v < 0, and so is V(w). 

Turning back from w to v, we first observe that both U(v) and V(v) are smooth in 
Vq < v < (and finite at v — 0). That is, no new singularity appears in the range v < 0. We 
still need to check the effect of the perturbation on the features of the singularity at v — 0. 
From Eq. (13) it is obvious that both V v and V vv are unaffected at the leading order in w, 
so the asymptotic behavior at v — is still correctly described by Eqs. (8,9). It then follows 
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that PP Riemann components [with respect to the tetrad (11)] diverge just as in the basic 
solution, e.g. 

e <* e P e l J E> _ O p (U-V) D „ U,|l/"-2 

(However, all curvature scalars vanish, as we are still dealing with a plane- wave solution.) 
The hypersurface v = thus remains a non-scalar PP curvature singularity. Since both U 
and V are finite at v = 0, the singularity is weak. We conclude that the non-scalar null 
weak PP curvature singularity of the basic solution is stable to small (but generic) ingoing 
plane-wave perturbations of the type considered here. 



IV. STABILITY TO OUTGOING PLANE-SYMMETRIC PERTURBATIONS 

Next, we check stability with respect to outgoing perturbations; namely, we shall now 
assume that non-trivial initial data are present also on the ingoing initial hypersurface v — Vq. 
The geometry is still plane-symmetric, but is no longer a plane-wave; Rather, it is described 
by the (LP) CPW solution. The line element is 



ds z = -e~ M dudv + e~ l 



y dx 2 + e- y dy 2 , (15) 



where the functions M, U, V generally depend on both v and u. This line element and the 
corresponding field equations have been discussed by several authors pT|p^| ; Here we shall 
briefly present the field equations in a form suitable for our analysis (basically the same form 
as in Ref. pi). The vacuum Einstein equations include five non-trivial equations, 



U uv - U U U V = (16a) 

M, uv - \{V U V V - U U U V ) = (16c) 

and 

2 U vv - U/ - V/ + 2M iV U v = (17a) 

2 U uu - U tU 2 - V u 2 + 2M, U U U = . (176) 

The evolution equations (16a-c) form a closed system of hyperbolic equations for the three 
unknowns U, V, and M. In the characteristic initial-value setup, the values of these three 
functions are to be specified on the two characteristic initial null hypersurfaces, u = uq and 
v — v o (see Fig. 2). We shall denote the initial values of U, V, and M on u = u by v U(v), 
v V(v), and v M(v), correspondingly, and those on the other null hypersurface, v — Vq , by 
u U(u), u V(u), and u M(u). We shall fix the gauge by demanding 

u M(u) = , v M(v) = . (18) 

Any solution of the evolution equations (16a-c) will satisfy the equations (17a,b), provided 
that the latter equations are satisfied on the two initial null hypersurfaces. Thus, in this 
gauge the Einstein equations are reduced to a set of three evolution equations (16a-c), 
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supplemented by the demand that the six initial functions v V(v), v U(v), v M(v) and u V{u), 
u U(u), u M(u) will satisfy Eq. (18) and the two ordinary differential equations, 

2 V U VV - V U/ - V V/ = (19a) 

2 U U UU - U U U 2 - U V U 2 = 0. (196) 

Correspondingly, a (LP) CPW spacetime is characterized by two arbitrary initial functions, 
v U(v) and u U(u); The two other nontrivial initial functions, v V(v) and u V(u), will be 
determined from Eqs. (19a,b). 




FIG. 2. The characteristic initial-value setup for the CPW solutions constructed in Sec. IV. 

In the class of CPW solutions, the ingoing plane-wave solutions form a subclass of mea- 
sure zero, characterized by a trivial initial function u U(u) = const. One then immedi- 
ately obtains from Eq.(19b) and Eqs.(16a-c) that U{u,v) = v U(v), V(u,v) = v V(v), and 
M{u,v) = 0, and Eqs. (1,2) are recovered. The plane-wave solutions considered in sec- 
tions II and III correspond to v U(v) = Uo(v) and v U{v) = U$(v) + AU(v), respectively 
(and constant U U). Here we shall consider nontrivial initial functions on both initial null 
hypersurfaces. On the outgoing hypersurface u = u , we take v U(v) = U (v) + AU(v) as 
in section On the ingoing hypersurface v = Vo, the initial function u U(u) can be any 
smooth function of u in the range uq < u < u\ (for some u\ > uq) which satisfies the 
following two obvious requirements: 

(i) 2 u U jUU — u U iU 2 > in «o < h < u\\ This will ensure the existence and smoothness of the 
function u V(u) [determined from Eq. (19b)] on the initial hypersurface v = Vq, throughout 
the relevant range Uq < u < u±. 

(ii) u U{u = uq) = v U(v = Vq), which is dictated by continuity of U(u,v) at the intersection 
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point (u ,v ). 

Our goal is now to analyze the evolution of geometry inside the domain of dependence, 

{u < U < U\ , Vq < V < 0}. 

We start by analyzing U(u,v). Defining 

W = e~ u , (20) 

Eq. (16a) is reduced to W >uv = 0. The solution, subject to the above initial data, is 

W(u, v) = v W (v) + u W (u) - u W (u ) , (21) 

where U > V W = exp(— u > v U) and 

°W = u W{u = u ) = v W(v = v ) . 

From the smoothness of v U(w) and u U(u) it then follows that W(u,w) is C°° in the entire 
domain u < u < u\ , v < v < 0, and hence W(u, v) is C°° in u < u < U\ , v o < v < 
(but not at v = 0; see below). Since U = — ln(W), the only possible singularity of U at 
v < is at the points where W vanishes. Obviously, 

W{u = u Q ,v) = v W(v) = exp [- v U(v)\ 

is strictly positive at Vq < v < 0. | By continuity of W(u, w), there exists a constant 112, 
Uo < U2 < ux, such that W > in the entire domain S = {u < u < u<i , Vq < v < 0} (see 
Fig. 2). Consequently, U(u,w) is smooth in S 1 , and U(u,v) is continuous in S 1 . In addition, 
U (u, v) is C°° in the domain S* - = {u < u < u 2 , fo < f < 0}. 

Before we analyze the evolution of V and M , we denote a remarkable feature of the 
evolution equations (16a-c): The form of these equations remain unchanged if we replace 
the independent variables u and v by new ones, u(u) and v(v) (with the unknowns U, V, 
and M unchanged); The only modification required is to replace du by du and dv by dv. 
We use this freedom to transform from v to w. Consider first Eq. (16b) (with v replaced 
by w). Since U(u,w) is known from Eqs. (20,21), this is a hyperbolic linear equation for 
V(u,w). Both its coefficients (i.e. U )U and U )W ) and its initial data [i.e. u V{u) and v V{w)] 
are C°° functions of u and/or w. Therefore, this equation has a unique smooth solution 
V(u, w) throughout S. Finally, consider Eq. (16c) for M(u, w). The second term in that 
equation, which includes derivatives of V and U with respect to u and w, is smooth in S, and 
the initial data at both initial hypersurfaces are M = 0. Consequently, M(u, w) is smooth 
throughout S. 

We conclude that all three functions U(u,w), V(u,w), and M(u,w) are C°° throughout 
S (even at w = 0). Returning now from w to the original independent variable v, we find 
that U(u,v), V(u,v) and M(u,v) are continuous throughout S and, moreover, are smooth 
in S~ . However, these functions will generically fail to be smooth at v — 0. To analyse 
this lack of smoothness, we use Eq. (7) to evaluate the v-derivatives of U, V and M. The 
maximal possible divergence rate of the first-order w-derivatives is l^l 1 /™ -1 (in fact, for U )V 



3 Note that V U is finite at v = 0, despite the divergence of V U <V there. 
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this divergence rate is never realized, because U }W always vanishes along the line v = 0, but 
this is not important for our discussion). In addition, as long as V w is nonzero at v = 
(which is indeed the case, as we shall immediately show), V vv is dominated by the last term 
in the right-hand side of Eq. (7): 



V 



vv — 



n- 2 (n-l)V w \v\ 1/n ~ 2 . (22) 



At u — «0j y.w is nothing but V V W , which is given in Eq. (13). Q Therefore, at u = uq we 
have 

V w {v = 0) = 2a(n- 2) / ( u = u ). 
From the smoothness of V(u, w) it then follows that 

V w {v = 0) ^ (23) 



in a neighborhood of u = Uq. I n this neighborhood, V vv thus diverges like \v 
The Riemann component 



l/n-2 



= (-v; w + £/^„ - M„y„) (24) 

is dominated at v = by V^„: 

= -^e^^)F w - \ [n- 2 (n - l)e^-^^] l^ 1 /- 2 (25) 

(see footnotes 1,2 above). Along an outgoing null geodesic of constant x, y, u, it is convenient 
to use the parallelly-propagated tetrad 

to evaluate the parallelly-propagated Riemann component 

<%4 e1e s x R aPlS = 2eV~ v+2M ) R vxvx - [ n - 2 (n - l)e 2A %] l^l 1 /-^ . ( 26 ) 

From Eq. (23) and the continuity of M it is obvious that at v = the term in brackets is 
nonvanishing in a neighborhood of u = uq. In this neighborhood, the parallelly-propagated 
Riemann component (26) diverges like \v J 1 /"- -2 . 

We find that in the perturbed spacetime, too, the hypersurface v — is a null curvature 
singularity. This singularity is weak, because U, V, and M are all finite at v — 0. We shall 
now show that this singularity is scalar-curvature. To that end, we shall calculate the scalar 



4 Despite the difference in the context, v V(v) is essentially the same as V(v) of section III, and, 
more specifically, their ^-derivatives (and hence also ^-derivatives) are identical. To verify this, 
note that V V )V is uniquely determined, through Eq. (19a), from the initial function v U{v) - in the 
same way that in the context of section |TJ V v is uniquely determined from U(v) through Eq. (2). 
Since v U(v) is identical to the function U{v) in Eq. (12), it follows that V V V is identical to E„ of 
Eq. (13). 
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K = R af3lS R a ^ S - Although the full expression for K is fairly complicated, a straightforward 
calculation shows that near v = it is dominated by 



K = 32 



2(M+U—V) p d , 2(M+U+V)r, d 

° 1 L vxvx A L uxux i ° A L vyvy A L uyuy 



This is the only part of K that includes second-order ^-derivatives, and it diverges like v x l n ~ 2 
(see below), whereas all other parts diverge like t> 2 / n_2 or slower. One can easily verify that 
the two terms in the brackets are equal, and that R uxux is of the same form as R vxvx [Eq. 
(24)], except that v is replaced by u. Expressing the dominant part of R vxvx in terms of 
Vvv, using Eq. (25), we obtain 

K = HV VV 

where 

H = 16 e 2M (V uu - U U V U + M iU V u ) . 

In appendix A we show that for a generic choice of initial functions u U(u) and AU(v), H is 
nonvanishing at v = in a neighborhood of u = u . In view of Eqs. (22,23), the scalar K 
diverges there like v 1 / n ~ 2 . The hypersurface v = is thus a scalar-curvature singularity. 

Let us summarize the main results of this section. The perturbed spacetime (described 
by a CPW solution) has the following features: 

i) No new singularity (spacelike or whatsoever) forms in S~ (i.e. before v = 0); 

ii) The hypersurface v = remains a null, weak, curvature singularity; 

hi) The divergence rate of the most divergent PP Riemann components is unchanged: v l l n ~ 2 . 
iv) For a generic outgoing perturbation, however, the singularity becomes scalar-curvature. 



V. DISCUSSION 

We have shown that within the framework of (linearly-polarized) plane-symmetric solu- 
tions the null weak singularity (4,5) is stable to both ingoing and outgoing perturbations: 
Both types of perturbations preserve the null weak character of the singularity. However, 
the outgoing perturbations generically convert the original non-scalar PP curvature singu- 
larity into a scalar-curvature singularity. This behavior is compatible with what we know 



about the CH singularity in black holes: In the mass- inflation model ||12|| , outgoing radi- 
ation converts the non-scalar PP singularity of the charged Vaidya solution ]T6]f27[ ] into a 
scalar-curvature singularity. Note also that the CH singularity in a generic spinning black 
hole is scalar-curvature. 0;i3 

The local stability of null weak singularities, demonstrated here (and also in Ref. |[25|| ), 
provides a strong support to the new picture of the CH singularity, which was obtained pri- 



marily from the perturbative approach [|K],|Tl[]. The present model of null weak singularities 
in CPW spacetimes may also serve as a useful toy-model for analysing various features of 
null weak singularities. 

Our analysis rules out the possibility that the introduction of outgoing perturbations 
will transform the entire null singularity into a spacelike one. It is still possible that, as the 
result of the outgoing perturbations, the null singularity will terminate at some u = u s > uq, 
where it intersects a spacelike singularity. This would be consistent with our construction, 
because, in the analysis in section 0, the demonstration of the regularity of U{u, w) is 
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restricted to the region S, i.e. u < u%. It thus may be possible that a spacelike singularity 
will be present at u > u-i- In fact, a spacelike singularity will positively form if W vanishes 
at some u < u < ui, as U = — hi(W) will diverge there. (Our analysis only guarantees that 
this cannot happen in the neighborhood of the initial hypersurface u = Uq. Note also that 
if the outgoing perturbation is bounded in a suitable way, then such a divergence will be 
excluded in the entire range uq < u < u\.) The line W(u, v) = (when exists) is known to 
be the locus of a spacelike, asymptotically Kasner-like, curvature singularity Note that 
this situation of a null singularity becoming spacelike at a certain point also occurs in the 
model of a spherically-symmetric charged black hole perturbed by a self-gravitating scalar 
field ||15|| . (At present, however, it is unclear whether such a situation also occurs in vacuum 
spinning black holes.) 

The solutions constructed here are all of the asymptotic form V oc v 1//n for some constant 



n > 2. This is also the situation in the analysis by OF [25[ - except that in the latter, unlike 



here, n was an integer. We shall refer to this type of asymptotic behavior as the power-law 
asymptotic behavior. The asymptotic behavior at the CH singularity of realistic black holes 
(as emerges from the perturbation analyses) is somewhat different, as we now explain: 

* For axially- symmetric perturbations of a Kerr BH [Kj (and also for a perturbed RN BH 
|T3|]), the asymptotic behavior is V oc (lnt>) _n , where n is the integer characterizing the 
power-law tails. We shall refer to it as the logarithmic asymptotic behavior. 

* For nonaxially- symmetric perturbation modes of a Kerr BH, the asymptotic behavior is 
V oc (ln|f|)~™ cos(mu; ln|w|), where u is a constant and m is the magnetic number of 
the mode in question. [|l(| We shall refer to it as the oscillatory-logarithmic asymptotic 
behavior. 

[Here v is the afline parameter along an outgoing null geodesic, with v = at the CH 
singularity, and V stands for a typical metric function (in a suitable gauge).] One would 
certainly like to extend the present analysis to these more realistic types of asymptotic 
behavior. It seems that the generalization to the logarithmic asymptotic behavior will be 
almost straightforward if one replaces w = \v \ l ^ n of the present analysis by w = (ln\v |) _n . 
The generalization to the oscillatory-logarithmic case is less straightforward, because of its 
non-monotonic nature. [In particular, when analysing Eq. (16b), it will not be possible to 
use w = (ln|f|) _n cos(u> ln|t>|) as a coordinate instead of v.] Still, due to the simplicity 
of the equations describing CPW spacetimes, hopefully the generalization to this type of 



asymptotic behavior will not be too difficult. Note that the approach used by OF [25 
seems to be inapplicable even for the logarithmic asymptotic behavior, because it is based 
on analyticity: Whereas for w = (— v) 1 ^ (with integer n) v is an analytic function of w, for 
w = (In H) - ™ v (w) fails to be analytic at v = 0. For the same reason, in Ref. |Z5| the power 
index n had to be an integer, whereas here n can be any real number > 2. In this respect, 
the present approach has an advantage over that of OF (however, as we mentioned in the 
introduction, in other respects the analysis by OF yields much more powerful results). 

The present construction was restricted to n > 2. It is straightforward, however, to 
extend it to any positive n ^ 1,2. The features of the resultant singularity at v — will 
significantly depend on the value of n. This issue deserves further investigation. It will also 
be interesting to generalize the present analysis to the logarithmic and oscillatory-logarithmic 
cases described above, and to the arbitrarily-polarized case. 
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APPENDIX A 



Let us define 

H (v) = (V uu -U u V u + M, u V u ) u=uo . 
It is sufficient to show that 

H (v = 0)^0. (Al) 

Then, H(u — u , v — 0) = 16e 2M H (v = 0) is nonvanishing too; and from continuity 
H(v = 0) is nonzero in some neighborhood u > u . 

In order to calculate Hq(v) we must evolve the ^-derivatives of U,V, and M along the 
line u = Uo, from v — vq (where these derivatives are obtained directly from the li-dependent 
initial functions), and up to v — 0. To that end we shall use the evolution equations (16a-c), 
which may be viewed as ordinary differential equations for the ^-derivatives of the metric 
functions. Let us define 

U u (v) = U tU {u = u ,v) , V u (v) = V u (u = u ,v) , M u (v) = M, a (u = u ,v) 

and 

V uu (v) = V uu (u = u ,v) , 

so 

H = Vuu - U U V U + M U V U . (A2) 
From Eqs. (20,21) it follows that 

U u = exp [ v U(v) -° U] u U u (u ) , 

where °U = v U(v — vo) = u U{u = u ). For our purpose, it is sufficient to recall that the 
functional dependence of U u on the initial functions is 

U u = U u [ v U; u U u (u )} . (A3) 

Here and below, the semicolon distinguishes between the functions of v (i.e. the initial 
functions at u — uq) at the left, and the parameters (obtained by evaluating the ^-dependent 
initial functions, and/or their derivatives, at u — uq) at the right. (For brevity, we omit 
the obvious dependence on v from this list of dependencies.) Next, applying Eq. (16b) to 
u = u , we obtain 

V = -(U V V + V U V ) 

V u (v) is to be determined from this ordinary differential equation, together with the initial 
condition V u (v = vo) = u V^ u (u ). It is straightforward to solve this linear equation explicitly 
(recall that U u , V U, and V V V are known), but again, for our purpose it is sufficient to recall 
that the functional dependence of V u on the initial functions is 

V u = V u [U u , V V V , V U V ; u V u (u )} ee V u [ V U ; u V u (u ), u U, u (u )} . (A4) 

[Here and below, we take into consideration Eq. (19a), which allows us to express V V V in 
terms of the function v U(v)\. In a similar manner, applying Eq. (16c) to u — Uq, we obtain 

M = -( V V V — V U U ) 

lvl u,v — 2\ v ,v v u u ,v u u) ■ 
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M u can thus be obtained by a direct integration, meaning that its functional dependence is 

M u = M u [V U ,U U , V V V , V U V ] ee M u u \/> ), u U u (u )} (A5) 

(recall that M u vanishes at v — vq). In summary, the three functions U u , V u and M u are 
fully determined from the function V U and the two parameters u V u (u ), u U jU (u ). 
Let us now analyse V uu . Differentiating Eq. (16b) we find 



U=UQ 



(A6) 



Using Eq. (17b), we obtain a differential equation of the form 

V uu ,v = \ V U V V UU + l -G{v) 
where 

G{v) ee [{U u 2 + V u 2 ) /2 - M U U U ] V V V + U U V U , V + U U , V V U . 
Note the functional dependence of G(v), 

G(v) = G[V U ,U U ,M U , V V V ] = G[ V U; u V u (u ), u U u (u )} . 

Equation (A6) is a linear ordinary differential equation for V uu (v), whose general solution is 

V uu (v) = exp ("U/2) \C+ l -"j G{v') exp [- v U(v')/2] dv' 

V vo 

where C is an integration constant. Recalling the initial condition, 

V uu (v = v ) = u V uu (u = u ) , 



we find 
where 

namely 



V uu (v) = Z(v) + exp [( V U - °U)/2] u V uu (u ) , 
Z(v) ee - exp( v U/2) J G(v') exp [- v U(v')/2] dv' , 



(A7) 



Z(v) = Z[ V U,G] = Z [ V U ; u V u (u ), u U u (u )} . 
Collecting Eqs. (A3,A4,A5,A7) and substituting in Eq. (A2), we can reexpress H as 



where 



H (v) = Q(v) + exp [( V U - °U)/2] u V uu (u ) , 
Q(v) ee Z(v) - U U V U + M U V U 



has the functional dependence 



Q = Q[ v U; u V u (u ), u U u (u )] 
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At v — 0, we find 



where 



H (v = 0) = g + exp {[ v U(v = 0) - °C/] /2) u V uu (u ) 



q = Q(v = 0) = g [^5 *%(«(>)> U U u ( u o)} ■ 



(A8) 



We still need to express the parameters u V u (uo) and u V uu (uo) in terms of the initial 
function u U(u), through Eq. (19b). The latter implies 



V,u — i ^2 U tUU U tU ^ , V^ un — i ( U tUUU U tU U tUU ) (2 U tUU U tU ^ 
so we can reexpress the functional dependence of go as 

q = qo[ v U; u U u (u ), u U uu (u )} . 
We can now write Eq. (A8) in the form 

H (v = 0) = ho + h u U, uuu (u ) , (A9) 

where 

1/2 



-1/2 



/i = go ± exp 
and 



[-u(v = 0) - °c/]/2] (- u c/, u M u c/, uu M) (2 u f/, uu ( Mo ) - tt f/,uW 2 )' 



-1/2 



/ii = ±exp = 0) - °U]/2\ [2 u U uu (u ) - u U u (u ) 2 

are two numbers with the same type of functional dependence: 

V = h 0>1 [ v U; u U u (u ), u U uu (u )} . 

Obviously, h\ 7^ 0. Therefore, in view of Eq. (A9), for any choice of initial function v U(v) 
and parameters u U tU (u ) , u U^ uu {uq) there only exists a single value of u U tUUU {uo) for which 
H (v = 0) = 0. Thus, for a generic choice of the initial function u U(u), the inequality (Al) 
is satisfied. 
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